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Lorentz-Dirac equations are currently accepted as the equations of mo- 
tion for charged particles in classical electrodynamics. These equations were 
derived by Dirac in a classical paper in 1938 the same equations were also 
obtained by Rohrlich on the basis of an action principle in his well-known 
book . In the case of two charges, the equations are: 

mi< = (ei/c) F^-vl + (e^/c) F^tvl + {H/^c') " (V^') a\a\v^) (1) 

where we use Rohrlich's notation. The equation for charge 62 is the same as 
(1), but with indexes 1 and 2 interchanged. The first term on the right-hand 
side (RHS) in (1), is Lorentz's force due to the external field. The second 
term on the RHS, links charge 62 charge 62 by means of a purely retarded 
interaction, and represents the only mutual interaction of the charged par- 
ticles in Lorentz-Dirac equations. Finally, the third term is the well known 
radiation reaction term, constituted by two parts: "Schott term", involving 
the third derivative of the charge position, and "Larmor term" , involving the 
square of the acceleration. In the case of only one charge in an external field, 
the second term on the RHS of (1) disappears. 

Equations in (1) are the final achievement of Dirac's conception, accord- 
ing to which it is possible to formulate a consistent set of classical equations 
of motion for point charges. Dirac reached these equations starting from 
Maxwell equations and the principles of relativity, masking the divergences 
associated with the singularities of the field at the position of the charges, 
by means of a mass renormalization procedure. Extensive theoretical work, 
carried out especially by Rohrlich [Q], led to consider that Lorentz-Dirac equa- 
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tions coherently describe the one-charge case in an external field. However, 
the situation does not remain the same when several charges are involved. In 
fact, a few years after Dirac derived his equations, Eliezer applied them 
to the case of two equal-mass particles with charges of equal magnitude but 
opposite signs, moving in a straight line and colliding head-on. This led 
to a controversy that still holds true [Q. According to Eliezer, the charges 
stop before colliding, then turn back and move away from each other, with 
an acceleration that is always different from zero, and with a velocity that 
tends toward the speed of light. This behavior is contrary to what common 
understanding would suggest for two charges of opposite signs, and is also 
contradictory to the law of energy conservation. Eliezer's conclusions were 
questioned by Clavier |p, Plass 0, and Rohrlich [0]. However, more recent 
detailed numerical calculations by Baylis and Huschilt [0, and Kasher [^], 
came to support Eliezer's conclusions. 

The pathological behavior of Lorentz-Dirac equations in Eliezer head-on 
collision has received, in general, little attention in textbooks (one exception 
is Parrott's book ). Several reasons may help to explain this lack of at- 
tention. In spite of the one-dimensional character of Eliezer's collision, the 
retarded nature of the interaction between the charges gives rise to mathe- 
matical complications. In particular, the numerical treatment of [0 & are 
rather elaborated, and questions about the convergence of the iterative pro- 
cess cannot be easily answered. This problem is related, in turn, to the issue 
of the existence and the uniqueness of the solution of Lorentz-Dirac equa- 
tions, which still remains an open question in the case of several charges. 
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Another source of troubles comes from the fact that in Ehezer's colhsion, 
the kinetic energy of the charges, the total rate of the radiation emitted by 
them, as well as the energy stored in the field of the charges, are all time- 
dependent. This makes it very difficult to analyze the energy conservation 
law at any given moment. Although the kinetic energy and the total rate 
of the radiation are well defined concepts, the same does not occur with the 
concept of "energy stored in the field of the charges" , which makes for even 
further troubles. We also want to point out that the catastrophic behav- 
ior of Eliezer's solutions has been obtained merely as a mathematical result 
of Lorentz-Dirac equations, but without identifying any physical cause- an 
identification that would be convenient for a deeper understanding. 

From our point of view, Lorentz-Dirac equations cannot be correct, be- 
cause they do not take into account all the radiation emitted by the charges. 
In fact, in the case of two charges, radiation must be described by: a Larmor 
term for each charge, terms that are indeed present in (1); plus a term that 
considers the interference effect between the fields of both charges, which 
is not present in (1). The second term on the RHS of (1) is the only one 
that considers both particles, but this term has very little to do with the 
interference radiation, which must involve the product of the acceleration of 
two charges, and this is clearly not the case. So, in order to put this con- 
jecture in a quantitative way, it would be ideal to find one exact solution of 
the Lorentz-Dirac equations for two charges where the total rate of radiation 
(including the interference term) can be calculated in a precise way. It is 
the purpose of this note to present one such solution, which is the first exact 
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solution for the Lorentz-Dirac equations for more than one charge present in 
the hterature related to the subject . As it will be shown below, Lorentz- 
Dirac equations, with appropriate external fields to be determined, allow for 
the motion of two charges moving in a plane at opposite ends of a diameter 
revolving at constant angular velocity in a fixed circular orbit of radius 
a. This solution can be considered as a sort of generalization of the exact 
solution found by Sokolov and Kolesnikova [|I^ for one charge in a circular 
orbit. 

Before discussing the algebraic treatment of our solution, it may be con- 
venient to make some general comments on Lorentz-Dirac equations, since 
treatment of this subject is almost absent in current literature. Besides, they 
may help those readers not very familiar with radiation reaction effects, and, 
on the other hand, they will allow us to set the framework of our discussion. 

It may seem rather natural to study the validity of Lorentz-Dirac equa- 
tions by means of experiments, instead of resorting to theoretical disquisition 
as we are doing here. However, due to other effects that appear in experi- 
ments that make it almost impossible to reach a clear conclusion, the troubles 
with Lorentz-Dirac equations cannot be clarified with the help of practical 
measurements. The motion of an electron in a storage ring is a good exam- 
ple for illustrating this situation, since in these machines radiation reaction 
plays an important role. The starting point for the description of the elec- 
tron motion in a storage ring is Lorentz-Dirac equations for one electron 



11, 0]. Now, in a machine of a few Gev, the magnetic external force is 



5 



much greater than the Larmor term, which in turn is much greater than the 
Schott term. These differences make it difficult to test Lorentz-Dirac equa- 
tions; in addition, there are several other complications. For example, in the 
vacuum chamber of the machine there is not only one electron, but many 
of them interacting with each other. In particular, these interactions lead 
to a loss of electrons from the beam. Furthermore, the particle trajectory 
is influenced by the interaction of the electron with the wall of the vacuum 
chamber, as well as by the inhomogeneity of the magnetic field. But the most 
striking phenomenon that rules out the use of Lorentz-Dirac equations for 
an accurate description of the motion are well-established quantum effects 



on electron motion Wu 12 



An experimental set-up of Eliezer's problem considering a collision be- 
tween an electron and a positron, will pose difficulties of the same kind as 
those observed in the case of a storage ring. So, even though Eliezer's prob- 
lem has a clear meaning in a classical context, where position and momentum 
have precise values at any time, quantum effects are unavoidable in the ex- 
periment and, therefore, Lorentz-Dirac equations do not apply. 

Why then should we insist on looking for exact solutions of Lorentz-Dirac 
equations if they are of little help for an accurate practical description of the 
motion of the charge? In our opinion, the study of these equations in their 
own context, regardless of experiments or practical applications, is interesting 
because their coherence is deeply related to the mass renormalization prob- 
lem, which still remains somewhat obscure in classical as well as quantum 
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theory. It was precisely in Lorentz-Dirac equations where, for the first time, 
Dirac introduced the concept of mass renormahzation, in order to deal with 
the divergent self-energy of a point charge. Thus, it would not be surprising 
to find that the violation of the energy conservation law in these equations 
is related to the mass renormalization problem. Now, if the Lorcntz-Dirac 
equations arc incorrect, as we clearly show below, the natural question that 
arises is: What are the correct equations of motion for point charges? This 
question assumes that the divergence due to the point nature of the charges 
can be consistently renormalized with the energy conservation law; but, is 
it possible that this program could not be carried out? In order to answer 
these and other related questions, Lorentz-Dirac equations must be studied 
in their own classical context. The analysis of exact solutions becomes a 
powerful tool for this purpose. 

Let us now describe our solution. We assume that the motion occurs on 
the X — Y plane, clockwise in a circle of radius a centered at the origin. 
The two identical particles of charges e < and mass m, are rotating with 
constant angular velocity ui, at the end of a diameter. At time t, let charge 
1 be the one located ai x — —a cos cut, y — a sin cut. Then, the retarded 
position of charge 2, which determines the field over charge 1 at time t, is 
defined by the retarded time t', time that also defines the retarded position of 
charge 1 respect to charge 2 at time t. This property is due to the symmetry 
of the motion. Now, let 2ip = cot — cot' be the angle between the diame- 
ters of the actual and the retarded positions. Then, we have (/? = /3cos(/?, 
where (5 — acu/c. In particular, we obtain dp/d(fi > 0. Thus, (3 is a strictly 
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monotonous increasing function of ip. Therefore, since there exists a one to 
one correspondence between /? and </?, we will use either of them to represent 
the velocity of the charges. Moreover, as /3 < 1, then, (/? < (/7cri = 0.739, and 
since /3 > 0, the angle </? is in the interval < </? < </7cri- 

In order to account for the retarded interaction between the charges 
(the second term on the RHS of (1)), we introduce the radial and tan- 
gential unitary vectors to the orbit at the location of charge 1, which are 
f = — i cos + J sin and t — i sina;^ + j cosa;i, respectively. Then, the 
radial component E^. — f ■ E and the tangential component Et — t ■ E oi 
the retarded electric field of charge 2 at the location of charge 1 at time t, 
are obtained from Lienard-Wiechert's well-known formulas for the field of a 
point charge. We get: 

(4a^/ e)Ej. — (cos (/? + (/p sin ip)~^ | {^ip sin 2ip -\- cos^ {l -\- </?^ cos 2ip sec^ (p^ 
— 2(/3^ (1 -|- tan cos 2(^1 , (2) 

(Aa^/e^Et — (cos(p + ipsmip)~^ {(ipcos2ip — sinipcosip) 

(l + </?^ cos 2(/? sec^ ip^ + 2ip'^ (1 + iptanip) sin2(/?| . (3) 

For this type of motion, the magnetic field at the position of each charge 
has only one component along axis Z, given by: 

Bz — — sin ipEr — cos (pEf. (4) 

The motion under consideration cannot be obtained without external fields. 
We will assume here the existence of a time-independent external electric 
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field tangent to the orbit, and with a fixed magnitude around the orbit circle. 
In order to show that this electric field arise from Maxwell's equations, we 
consider the ideal sources: a charge density p(x, t) that vanishes identically 
everywhere, and a charge current vector proportional to time t given by 
J(x, t) = —At6{r — 6)0, where we use cylindrical coordinates (r, 0, 2;); 6 is 
the usual Dirac delta function, and (f) is the unit vector tangent to the circle 
r = fixed, pointing in the direction of increasing 0. Parameters A and b are 
positive, with radius b of the infinitely long solenoid smaller than orbit radius 
a. Then, it is easy to show that for these sources the fields: 

E(x,i) = {pio/2)A{r-e{r-b){r-byr)}4> 

B(x,t) = -noAt{l-e{r -b)}k (5) 

where d{r) represents the step function, are solutions of Maxwell's equations 
everywhere. The proof comes easier in cylindrical coordinates. This proof 
uses either standard rigor with step and delta functions, or the distribution 
theory. We need the fields (5) only in region r > 6, in which case they reduce 
to E = {ijLoAb^/2r) ^; B = 0. The tangential electric field can take any value 
at a given orbit circle. In addition to the external electric field, we will also 
consider an homogeneous time-independent external magnetic field pointing 
in the negative direction of axis Z. The fourth component of Lorentz- Dirac 
equation (1) is then identically satisfied by 2; = at any time. Equation 
(1) for the component along the X axis, x — —a cos cut, of charge 1, can be 
combined with the component along the Y axis, y = a sin cot, of the same 
charge, in order to write them in terms of radial and tangential component 
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as follows: 

t {eEt + e£;r* - (2eV3a2) P^i"^} = r {-mau;^7 - eEr + e^B, + e/3Sf*} , 

(6) 

where E^^^ is the tangential component of the external electrical field, Bf^^ 
is the external magnetic field along the negative axis Z, and 7 is the usual 
relativistic parameter given by (1 — It is easy to show that charge 2 

also satisfies equation (6). Prom this equation, it follows that must have the 
following value: 

eEr^(eya')m, (7) 

where function /(</?) is given by: 

/((^) = (2/3)(^3 cos (^(cos' - (^2)-2 _ {aye)Et. (8) 

A detailed analysis shows that f{(p) is a positive strictly monotonous increas- 
ing function of (f. For this reason, for a given radius a and a given velocity 
V = au, equation (7) determines a perfectly well-defined external field E^^^, 
which in turn imposes a restriction on parameters b and A in (5). 

The value of the external magnetic field obtained from (6), is the 
following: 



B 



ext 



{e/Aaro) I iip (cos^ ip — (f'^'j ^ +{ro/a)ip ^ (ia^ / E^ip cos ip 



+ (4:0^/ ej Er (cos (p + (psin (^) j | , (9) 

where tq = e'^/mc^. For given radius and velocity, the above formula de- 
termines Bf^^ in a unique way. For (p tending to zero in (9), Bf^^ tends to 
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infinity, so that the magnetic fields compensate the repulsive electric field 
force between the charges. 

The first component of (1), representing the energy conservation law, 
reduces to: 

2ev • E*^^* = (Ae^c/2>a') (3^ - {4.6^0/30") (3^g{ip), (10) 
where g{ip) is given by: 

g{ip) = (3aV2e) ip'^ cos^ (pEf (11) 

Except for one factor, equation (10) coincides with equation (7). The left 
side of equation (10) corresponds to the rate at which the external field sup- 
plies energy to the system composed of the two charges. Now, the kinetic 
energy of the charges remains unchanged in this motion. Besides, the energy 
stored in the field of the charges remains also unchanged, since the position 
of the charges at two different times cannot be distinguished. We empha- 
size that we need not pay attention to the precise definition of the concept 
of energy stored in the field of the charges, since by the symmetry of the 
motion, this energy -whichever it may be- is time-independent. Therefore, 
according to the energy conservation law, all the power supplied to the sys- 
tem of two charges by the external electric field must be radiated away. This 
means that the RHS of (10) corresponds to the total rate of radiation emitted 
by the two charges. But for known motion of the charges. Maxwell's equa- 
tions determine the field generated by the charges, which in turn determines 
the total rate of radiation for this motion. We carried out this somewhat 
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elaborated calculation in a separate paper |T3[ . The total rate of radiation 
involves a Larmor term for each charge, plus an interference term due to the 
field of both charges. The Larmor terms for the two charges are precisely 
the first term on the RHS of (10), so its second term must coincide with 
the interference radiation term. If we put the interference radiation term in 
[|l^ in the form — (4e^c/3a^) f3^I{(p); then, function /(v^) must coincide with 
function g{(f) in (11). We plotted these functions in Fig. 1, showing that 
for velocities of the charges not close to the speed of light, function g{if), 
which comes from Lorentz-Dirac equations, is qualitatively similar to 
However, these functions are quantitatively different for any velocity of the 
charges, and therefore Lorentz-Dirac equations violate the energy conserva- 
tion law. We remark that if we consider only one charge in circular motion. 



as in the Sokolov-Kolesnikova solution [10| , the contradiction with the en- 
ergy conservation law disappears, since in this case we have to take Et = 
in (11). 



Figure 1 

The above discussion suggests that the source of troubles with the energy 
conservation law in Lorentz-Dirac equations for more than one charge, comes 
from an inappropriate consideration of the interference radiation between the 
charges. So, it seems natural to solve this problem by incorporating the in- 
terference radiation energy in the equations of motion. There are motions 
for which this program may succeed, as seems to be the case for the motion 
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considered herein. Nevertheless, this procedure may not work in all cases. 
Charges in close collision may present troubles, because the trick of oversim- 
phfying the relevance of the fields in the neighborhood of the singularities of 
the charges by means of the mass renormalization procedure may critically 
disturb the exchange of energy momentum between the singularities. This 
disturbance may lead to a violation of the energy conservation law. 
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FIGURE CAPTIONS 



Figure 1: The solid line represents function /((/?) of the interference rate 
of radiation given by — (Ae^c/Sa^) P^I{(p). The dotted line represents the 
corresponding function g{<^) given by (11), for the interference rate according 
to Lorentz-Dirac equations. 
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